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ABSTRACT. Time-discrete dynamical systems on a finite state space have been used with great success to 
model natural and engineered systems such as biological networks, social networks, and engineered control 
systems. They have the advantage of being intuitive and models can be easily simulated on a computer in 
most cases; however, few analytical tools beyond simulation are available. The motivation for this paper is 
to develop such tools for the analysis of models in biology. In this paper we have identified a broad class 
of discrete dynamical systems with a finite phase space for which one can derive strong results about their 
long-term dynamics in terms of properties of their dependency graphs. We classify completely the limit 
cycles of scmilattice networks with strongly connected dependency graph and provide polynomial upper and 
lower bounds in the general case. 



1. Introduction and Background 

Time-discrete dynamical systems on a finite state space play an important role in several different contexts. 
Examples of such systems include Boolean networks, cellular automata, agent-based models, and finite state 
machines, to name a few. This modeling paradigm has been used with great success to model natural and 
engineered systems such as biological networks, social networks, and engineered control systems. It has the 
advantage of being intuitive and models can be easily simulated on a computer in most cases. A disadvantage 
of discrete models of this type is that few analytical tools beyond simulation are available. The motivation 
for this paper is to develop such tools for the analysis of models in biology, but the results are of independent 
mathematical interest. Some theoretical results have been proven for Boolean networks and are reviewed 
in [7J. In that paper the authors study conjunctive Boolean networks, that is, Boolean networks for which 
the future state of each node is computed using the Boolean AND operator. It is shown that the dynamics 
of such networks is controlled strongly by the topology of the network. The current paper shows that the 
results in [7J are valid much more broadly. We briefly describe the main results in [7J. 

A Boolean network / on n nodes be viewed as a time discrete dynamical system over the 

Boolean field F2: 

/ = (/!,...,/„) :FJ— >F5, 

where the coordinate functions fi are the Boolean functions assigned to the nodes of the network. We can 
associate two directed graphs to /. The dependency graph (or wiring diagram) of / has n vertices 

corresponding to the variables Xi, ■ ■ ■ ,x n of /. There is a directed edge i — > j if the function fj depends on 
X{ (i.e. there is an instance where changing the value of Xi changes the value of fj). That is, £>(/) encodes 
the variable dependencies in /. The dynamics of / is encoded by its phase space, denoted by S(f). It is 
the directed graph with vertex set F^ and a directed edge from u to v if /(u) = v. Thus, the graph T)(f) 
encodes part of the structure of / and the graph S(f) encodes its dynamics. The results in [7J relate these 
two graphs, deriving information about <S(/) from Q(f), in the case where each /, is a conjunction of the 
variables Xj for which there is an edge j — >• i in These networks were called conjunctive networks in 

It is easy to see that the graph S(f) has the following structure. Its connected components consist of a 
unique directed cycle, a limit cycle, with directed "trees" feeding into the nodes of the cycle, representing 
the transient states of the network, each of which eventually maps to a periodic point. If the graph D(f) is 
strongly connected, that is, there is a directed path from any vertex to any other vertex, then it is shown 
in [7j that there is a precise closed formula for the number of limit cycles of any given length. This formula 
depends on a numerical invariant of £>(/), its loop number. If S(f) is not strongly connected, then there is 
a sharp lower bound and an upper bound on the number of limit cycles, in terms of the loop numbers of the 
strongly connected components of S(f) and the antichains in the poset of strongly connected components. 
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These two bounds agree for the number of fixed points of /, so that one corollary is a closed formula for the 
number of fixed points of a general conjunctive network. 

The current paper shows that these results hold in much broader generality, namely for time discrete 
dynamical systems over any finite set X: 

f = (f 1 ,...,f n ):X n ~^X n , 

such that the fa are constructed from an operator A : X 2 — > X, which has the property that it endows X 
with the structure of a semilattice, that is, A is commutative, associative, and idempotent. We will call such 
functions fa semilattice functions. That is, any semilattice X gives rise to a dynamical system to which the 
formulas and bounds for limit cycles hold. It is shown furthermore that they hold precisely for the class of 
semilattice networks. It is important to mention that since semilattice networks are not linear (see Example 
16. ip , our results complement those for linear systems in [3J [5] . 

2. Semilattice networks 

Let X be a set with M elements and consider a dynamical system with n variables over X: 

f = (f 1 ,...J n ):X n —>X n . 

Definition 2.1. A function A : X 2 — > X is called a semilattice operator if it satisfies the following (with 
notation x A y := A(x, y)): 

x A y = ?/ A x , A is commutative 

x A (y A z) = (x A y) A z . A is associative 

x /\x = x , A is idempotent 

Note that semilattice operators induce the structure of a semilattice on X. Conversely, the meet or join 
operations on a semilattice are semilattice operators in the above sense. 

Example 2.2. Some examples of semilattice operators are: 

A = AND over [0, l] 2 
A = OR over [0, l] 2 
A = MIN over [0,m] 2 
A = MAX over [0,m] 2 

That is, semilattice operators are a generalization of the conjunctive and disjunctive operators used in [7]. 
The domain of a semilattice operator A : X 2 — > X can be extended to A : X k — > X with k > 1 by 
x\ A X2 A . . . A Xk — x\ A (x2 A (. . . A (xk-i A Xk))) if k > 1 and 
Ax = x A x = x 



We will call such an extended operator a semilattice function. 

Definition 2.3. A dynamical system / = (/i, . . . , /„) : X n — > X n is called a semilattice network if there 
exists a semilattice operator, A, such that fj = /\ is a semilattice function for all j = 1, . . . , n (where fj 
depends only on x^ , . . . , Xi k ). 

Let / : X n — > X n be a semilattice network, G = £>(/), and A the adjacency matrix of G. We will 
assume here and in the remainder of the paper that none of the coordinate functions of / are constant, that 
is, all vertices of G have positive in-degree. 

2.1. Structure of the Dependency Graph. Define the following relation on the vertices of G: a ~ b if 
and only if there is a directed path from a to b and a directed path from b to a. It is easy to check that ~ 
is an equivalence relation. Suppose there are t equivalence classes Vi, . . . , Vt- For each equivalence class Vj,, 
the subgraph Gi with the vertex set Vi is called a strongly connected component of G. The graph G is called 
strongly connected if G consists of a unique strongly connected component. A trivial strongly connected 
component is a graph on one vertex and no self loop. Since such components do not influence the cycle 
structure of the network, we assume that all strongly connected components are non-trivial. 
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Example 2.4. Our running example for this paper will be the semilattice network / : {0, 1, 2} 63 — > {0, 1, 2} 
with dependency graph given by Figure [Hand semilattice function A = MIN . 
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Figure 1 . Dependency graph of the semilattice network in Example 12.41 The labels have 
been omitted for simplicity. 



Let Gi be a strongly connected component, and let hi be the semilattice network with dependency graph 
D(hi) = Gi. Let h : X n — > X n be the semilattice network defined by h = (hi, . . . , h t ). That is, the 
dependency graph of h is the disjoint union of the strongly connected graphs G%, . . . , Gt, and h is obtained 
from g by deleting all edges between strongly connected components. 

Now define the following order relation on the strongly connected components G\, . . . , Gt of the dependency 
graph of the network /. 

(2.1) Gi < Gj if there is at least one edge from a vertex in G, to a vertex in Gj. 

In this way we obtain a partially ordered set V . In this paper, we relate the dynamics of / to the dynamics 
of its strongly connected components and the poset V . 

Example Q] (Cont.). The dependency graph of / has four strongly connected components, Gi, G2, G3, 
and G4 (bottom, left, right and top, resp.). The poset is given by G\ < G2, G3 -< G4. 




Figure 2. The strongly connected components of / (left) and their poset (right). 
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2.2. The Loop Number. 

Definition 2.5. The loop number of a strongly connected graph is the greatest common divisor of the 
lengths of its simple (no repeated vertices) directed cycles. The loop number of any directed graph G is the 
least common multiple of the loop numbers of its non-trivial strongly connected components. 

Example[T](Cont.). In our running example, the loop numbers are loop{Gi) — 1, loop(G 2 ) — 2, loop{G^) = 
3 and loop{Gi) = 1. 

3. Semilattice Networks with Strongly Connected Dependency Graph 

In this section we give an exact formula for the cycle structure of semilattice networks with strongly 
connected dependency graphs. In the next section we will also consider networks with general dependency 
graphs and give upper and lower bounds for the cycle structure. 

The formula for conjunctive networks was given and proven in [7]. It is not difficult to show that the 
proofs are also valid for general semilattice networks. 

Let / : X n — !► X n be a semilattice network with semilattice operator A. Assume that the dependency 
graph of / is strongly connected with loop number c. 

Lemma 3.1. The set of vertices of 2>(/) can be partitioned into c non-empty sets Wi, . . . , W c such that each 
edge of'D(f) is an edge from a vertex in Wi to a vertex in Wi+i for some i with 1 < i < c and W c +i = W±. 

Proof. For a proof of this fact see [TJ Lemma 3.4.1 (iii)] or [5J Lemma 4.7]. □ 

Let Si be the number of elements of Wi. Without loss of generality we assume for the rest of the section 
that Wi = {1, ... , si}, W 2 = {si + 1, . . . , si + s 2 },. . . , W c = {si + . . . + s n _i + 1, . . . , s 1 + . . . + s c = n}, 

Proposition 3.2. Let c be the loop number of f; then there exists ko such that for all k > k$ we have 

f ck (x) = ( yi,...,yi ,..., y c ,---,y c ) 

Si times s c times 

f ck+1 (x) = ( y c ,--j,yc , yi,..j,yi , y c -i, ■ ■■ , y c -i) 

si times S2 times s c times 

f ck+J (x) = ( y c +i-j, ■ ■ ■ , 2/c+i-j , • ■ ■ , y c , ■ • • , yc , yi, ■ ■■ , yi , ■ ■ -, y c -j,- ■ ■ , y c -j ) 

si times Sj times Sj + i times s c times 

r k+c (x) = r k (x) 

where y 3 = A, eWj x t 

Proof. The proof is analogous to [2j Theorem 4.10]. □ 

The following corollary states that the period of / can be obtained from the topology of its dependency 
graph. 

Corollary 3.3. The period of f is equal to the loop number of'D(f). 

The following corollary states that the long-term dynamics of a semilattice network can be reduced to the 
dynamics of a rotation over X c , with c as in Lemma |3. II 

Corollary 3.4. The cycle structure of f is equal to the cycle structure of R : X c — > X c where R(yi, ■ ■ ■ , y c ) = 
(y c ,yi, ■ ■ ■ ,2/c-l)- 

Proof. Let r : X c -> X n and $ : X n X c be defined by T(y 1 ,y 2 , ...,y c ) = ( yi, . .. , yi , y 2 , . .. , y 2 , . . . , y c ,---,yc) 

si times S2 times s c times 

and . . . , x„) = {x Sl ,x Sl j rS2 , . . . ,x„). The proof now follows from the equalities $ o / o T = t and 

r O T O $ = /| {periodic points of /} ■ □ 
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For any positive integers p, k that divide c, let A(p) be the set of periodic points of period p and let 
D(k) :=\jA(p). 
p\k 

Proposition 3.5. The cardinality of the set D(k) is \D(k) \ — M k . 

Proof. It follows from the fact that if k\c, then a rotation in M colors and c variables has M k colorings of 
periods that divide k. □ 

Corollary 3.6. If p is a prime number and p k divides c for some k > 1, then 

\A(p k )\ = MP" - MP k ~" . 

Proof. It is clear that |-D(1)| = M (there are M constant colorings). Now if p is prime and k > 1, then the 
proof follows from the fact that D{p k ) = D(p k ^ 1 ) {*)A(p k ), where (+J is the disjoint union. □ 

Next we prove Theorem 13. 71 which gives the exact number of periodic points of any possible length. 

Theorem 3.7. Let f be a semilattice network whose dependency graph is strongly connected and has loop 
number c. If c = 1, then C(f) = MC . If c > 1 and k = p^ 1 . . .p s r T is a divisor of c, then the number of 
periodic states of period k is \A(k)\ 

l l 



(3.i) \A(k)\ = J2 ■ ■ ■ J2 (~iy i+i2+ - +tr m p 



i 1= Q i r =0 

Proof. The statement for c = 1 is part of the previous corollary. Now suppose that c > 1. For 1 < j < r, let 
kj = pj 3 1 n[=i i=£jPT ■ Then D(k) = A(k) |+)(U^=i D(kj)), where 1+J is a disjoint union, in particular, 

r 

A(k) = D(k) \ [j D{k 3 ). 

The formula 13.11 follows from the inclusion-exclusion principle and the disjoint union above. □ 

\A(k)\ 

Corollary 3.8. If k divides c, then the number of cycles of length k in the phase space of f is C(f)k = — - — . 
Hence the cycle structure of f is 

c(/)=i; m c *. 

fc|c 

Example [1] (Cont.). In our running example we obtain: 

C(hx) = 3Ci, C(h 2 ) = 3Ci + 3C 2 , C(h 3 ) = 3Ci + 8C 3 , C(/i 4 ) = 3C X 

Remark 3.9. Notice that the cycle structure of / depends on the loop number and \X\ = M only. In 
particular, a semilattice network with loop number 1 on a strongly connected dependency graph only has as 
limit cycles the M fixed points (w, w, . . . ,w) where w € X , regardless of how many vertices its dependency 
graph has and how large n is. 

4. Networks with general dependency graph 

Let / : S — > S be a semilattice network with dependency graph £)(/). Let G%, . . . ,Gt be the strongly 
connected components of S)(/). Furthermore, suppose that none of the Gi is trivial. For 1 < i < t, let hi be 
the semilattice network that has Gi as its dependency graph and suppose that the loop number of hi is q. 
In particular, the loop number of / is c := lcm{ci, . . . , ct}. 

First we study the effect of deleting an edge in the dependency graph between two strongly connected 
components. Let G\ and G2 be two strongly connected components in 3D(/) and suppose G\ -< Gi. Let 
x — > y be a directed edge in £>(/) between a vertex x in G\ and a vertex y in Gi. Let 2D' be the graph 
£>(/) after deleting this edge, and let g be the semilattice network such that D (g) = D' . 

Theorem 4.1. Any cycle in the phase space of f is a cycle in the phase space of g. In particular C(f) < C(g) 
componentwise. 
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Proof. First notice that x appears in the expression f£ if and only if there is a path from x to y of length 
fc. Let C := {u, /(u), . . . , / m-1 (u)} be a cycle of length m in S(f). To show that C is a cycle in S(g), it is 
enough to show that, u x A uq 2 = Au<3 2 . Thus, the value of y is determined already by the value of y' and 
the edge (x,y) does not contribute anything new and hence C is a cycle in S(g). This is equivalent to show 
that u x A Uj,' = Uyi for all y' G G2 such that there is an edge from y' to y. 

Suppose the loop number of G\ (resp. G2) is a (resp. b). Now, any path from x (resp. y) to itself is of 
length pa (resp. qb) where q,p > T and T is large enough, see [2j Corollary 4.4]. Thus there is a path from 
x to y of length qa + 1 for any q > T. Also, there is a directed path from y to y' of length qb — 1 for any 
q > T. This implies the existence of a path from x to y' of length q(a + b) for all q > T, then x appears in 

•V ' 

Now u = / m (u) = / mfe (u), for all fc > 1. Choose q, k large enough such that q(a + b) = km > T. Then, 
IV = /™ fe ( u ) = u x A . . . and u x A ity = A u K A . . . = A . . . = u y >. Therefore u x A u G2 = Au G2 . □ 

Let h : S — > S be the semilattice network with the disjoint union of G\, . . . , Gt as its dependency graph. 
That is, h = (hi, . . . , h t ). For h, there are no edges between any two strongly connected components of 
the dependency graph of the network. Its cycle structure can be completely determined from the cycles 
structures of the hi alone. 

Theorem 4.2. Let C(hi) = ajjC be the cycle structure of hi. Then the cycle structure of h is 

= lli=i C(hi) (where C r C s := lcn ^. C lcm ^ r ' s ^ ) and the number of cycles of length m (where m\l ) in the 
phase space of h is 

(4.1) C(h) m = Y, ^-^"l l • 

lcm{j'i,...,j t }=m 

Proof. This follows from the fact that if u is a periodic point of hi of period fc,; and v is a periodic point of 
hj of period kj, then (u, v) is a periodic point of (hi, hj) of period lcm(fci, kj). □ 

Corollary 4.3. Let f and h be as above. The number of cycles of any length in the phase space of f is 
less than or equal to the number of cycles of that length in the phase space of h. That is C(f) < C(h) 
componentwise. In particular, the period of f is a divisor of the loop number of its dependency graph. 

In [7] , it was shown that the posct structure of V gives an algebraic way to combine the cycle structure of 
hi to obtain lower and upper bounds for the cycle structure of /, where / was a conjunctive Boolean network. 
It is not difficult to see that the proofs of these results still also hold for general semilattice networks if the 
corresponding semilattice operator has a "neutral" and an "absorbent" element (analogous to the identities 
1 A x = x and A x — 0). In order to state the theorem about lower and upper bounds of semilattice 
networks, we need the following definitions. 

Definition 4.4. Let A : X 2 ->Xbea semilattice operator. An element A G X is called a neutral element 
if A A x — Ax for all x G X. An element 9 G X is called an absorbent element if x A 6 = A6 for all x G X. 

Example 4.5. All the functions in Example 12.21 have a neutral and absorbent element; they are: 

A = 1 , 9 = 
A = , 9 = 1 
A = m , 9 = 
A = , 9 = m 

Remark 4.6. Since X is finite, every semilattice operator has an absorbent element (9 = A x exx). 

Remark 4.7. Any semilattice operator A on X can be extended to a set with a neutral element by defining 
A' : (X U {A}) 2 4lU {A} as x A' y = x A y if x, y G X and x A' A = A A' x = x otherwise. 

Let f,h,Gi,...,G t be as above and let Cj be the loop number of G^. Furthermore, let V be the poset of 
the strongly connected components. Let £1 be the set of all maximal antichains in P. For J C [f] = {1, . . . , t}, 
let xj := Ylj^jXj and let J := [t] \ J. In the remainder of the paper we assume that / : X n —> X n is a 
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semilattice network such that its semilattice operator, A, has idempotent neutral and absorbent elements, A 
and 6, resp. 

Definition 4.8. For any subset J C [t], let 

J- := {k : Gj r< Gk for some j e J}, J- := {fc : Gj y Gk for some j e J}, 

:= {k : Gj -< Gfc for some j € J}, and := {fc : Gj >- Gfc for some j e J}. 
A limit cycle C in the phase space of / is Jg (resp. J\) if the Gj component of C is 9 (resp. A) for all j G J. 

Denote (A', L) = j °' if ^ n L ^ and j = jfc y TV, J w = U M. 
X ' ' [1, ifXni = 0. 

Definition 4.9. The L- and {/-polynomial associated to V are defined as follows: 

c( Z1 ,..., 2t ) = n ** 



Zk 



iQNC[t] kei N uJ* 

JCMC[t] 



Example [T] (Cont.). For the poset in Figure[5]we obtain: 

C(zi,z 2 ,z 3 , Zi) = -2 + zi + z 2 z 3 + z A 

U(zi,z 2 , Z3, Zi) = 14 - 7z\ + iz\z 2 — 4z 2 + 3ziz 3 - ziz 2 z 3 + z 2 z 3 - 4z 3 + 4ziz 4 - 2z\z 2 za 
+ 3z 2 Z4 - 2ziz 3 z 4 + ziz 2 z 3 Zi - z 2 z 3 Z4 + 3z 3 z± - 7z 4 

Theorem 4.10. With the notation above we have the following coefficient-wise inequalities 

C(C(hi), . . .,C(h t )) < C(f) < U{C{h!), . . .,C(h t )). 

Here, the polynomials C and IA are evaluated using the "multiplication" described in Theorem an d 
coefficient-wise addition. 

Proof. The proof is analogous to the proofs of Theorems 6.2 and 7.4 in [7]. □ 

Note that the left and right sides of the inequalities (|4.10[) are polynomial functions in the variables C(hi), 
with integer coefficients. That is, the lower and upper bounds arc polynomial functions depending exclusively 
on measures of the network topology. 

Example [T] (Cont.). In our running example we obtain: 

C(C(hi),C(h2),C(ha),C(h A )) = 13d + 9C 2 + 24C 3 + 24C 6 
U{C{h x ),C{h 2 ),C{hz)X{hi)) = 20d + 24C 2 + 64C 3 + 96C 6 

and therefore 13Ci + 9C 2 + 24C 3 + 24C 6 < C(f) < 20Ci +24C 2 + 64C 3 + 96C 6 . It is important to mention that 
the phase space of / has 3 63 w 10 30 nodes so it is not feasible to obtain information about the cycle structure 
from exhaustive enumeration. Also, although the bounds agree on fixed points for Boolean semilattice 
networks [7], our example shows that they do not agree for general lattice networks. 

5. Characterization of Semilattice Networks 

In this section we characterize semilattice networks; in order to do this, it is enough to characterize 
semilattice operators. Since semilattice operators are semilattice operations, the number of semilattice 
operators over a set with m elements (up to permutation) is the number of semilattices with m elements. 
According to the next proposition, the number of semilattice operators over a set with m elements is the 
number of lattices of size m + 1. Although there is no closed formula for the number of lattices of a given 
size, algorithms for counting them have been developed [4]. 

Proposition 5.1. There is a one-to-one correspondence between semilattices with m elements and lattices 
with m + 1 elements. 
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Proof. If (X, A) is a semilattice with m elements, consider (X U {A}, A', V), by denning x A' y = x A y if 
x,y £ X and iA'A = AA'i = i otherwise; also, x V j/ = A{z :zAi = i, z Ay = y}. It follows that 
(X U {A}, A', V) is a lattice with m + 1 elements. On the other hand, if (Z, A, V) is a lattice with m + 1 
elements, let A = /\ Z; then it follows that (Z\{X}, A|^\{a} 2 ) is a semilattice with m elements. □ 

6. Infinite Semilattice Networks 

We present some results on infinite semilattice networks, both networks on an infinite set X and networks 
on an infinite Cartesian product of a set X. 

6.1. Semilattice networks on infinite sets. If X is a set with infinitely many elements (that is, M = oo), 
some of the theorems remain valid. Notice that if S C X is finite, there exists a finite set ZDS that is 
closed under A; that is, A\z ■ Z 2 — > Z is a semilattice function. 

Suppose that is strongly connected with loop number c, and let u = (ui, . . . , u n ) be a periodic point 
of /. Let Z D {ui, . . . ,u n } be a finite subset of X that it is closed under A; then we can consider u as a 
periodic point of f\z™- Therefore, the results on finite semilattice networks apply, and the period of u must 
divide c. That is, Corollarv l3.3l is valid for M = oo. Now, consider a divisor k of c and consider Z <Z X finite 
with at least m elements such that A\z? is a semilattice operator. Then, the number of periodic states of 
period k of f\z*> is at least |A(fc)| (sec Theorem 13 . 71 and notice that |A(fc)| is increasing with respect to M). 
Since limm-^ |^4(A;)| = oo, it follows that / has infinitely many periodic points and limit cycles of length k. 
Then Theorem 13.71 (and the corresponding corollary) holds for M = oo. 

If is not necessarily strongly connected, suppose A has a neutral and an absorbent element. Let 

hi, . . . , h n correspond to /|z™Then, by using the argument in the paragraph above, it follows that if at least 
a limit cycle of length k appears in C{C{h{), . . . , C(ht)), then / has infinitely many limit cycles of length k. 

6.2. Infinite-dimensional semilattice networks. If the dimension of / is infinite, that is, n = oo, then 
A needs to satisfy an additional condition to be properly defined: Every (possibly infinite) subset S C X 
has a largest lower bound. This means that there exists x £ X such that iAs = i for all s G S (in lattice 
terminology: x is a lower bound); and if there is another such x', then iAi' = x' (in lattice terminology: x 
is the largest lower bound). This additional condition allows for a function /, to have infinitely many inputs. 

Suppose that S)(/) is strongly connected with loop number c, and let u = (ui, . . .) be a periodic point 
of / of period d. Consider x,y <E Wi (see Lemma |3.1|) : then there exists kg such that for all k > ko 
there is a path of length ck from x to y and from y to x. Then, fy k = x A . . ., f x k = y A . . . for all 
k > k ; in particular, u y = fy k ° d (u) = u x A w and u x = f^ k ° d (u) = u y A v for some v, w £ X. Then, 
u x A u y = u x A u x A w = u x A w — u y , similarly u x A u y = u x ] therefore u x = u y . It follows that Corollary 
13.31 and Theorem 13.71 remain valid for n = oo (and M £ Z + U {oo}). 

If £>(/) is not necessarily strongly connected, suppose A has a neutral and an absorbent element. It is 
not difficult to show that Theorem 14. 101 remains valid for n = oo (and M G Z + U {oo}). 

Finally, we show with counterexamples that we cannot omit any of the properties of A in Definition 12.11 
That is, the formulas and bounds derived in this paper are valid exactly for semilattice networks. 

Example 6.1. Consider A = +, then A is commutative and associative, but not idenrpotcnt. Consider 
/ : — > F|, defined by f(xi, x%) = (xi + X2, xi + X2). It is not difficult to see that / has a unique limit cycle 
(a fixed point), (xi,X2) = (0,0); that is, C(f) = C . On the other hand, the loop number of its dependency 
graph is 1, so from Theorem 13.71 we would obtain 2C 1 ^ C(/). 

Example 6.2. Consider x Ay — x 2 y defined over F3. It is easy to show that A associative and idenrpotcnt, 
but not commutative. Consider / : F| — > F|, defined by /(xi,X2) = (x^X2,xi). It is not difficult to show 
that / has the limit cycle of length 2, {(1, 2), (2, 1)} (it also has 3 fixed points). On the other hand, the loop 
number of its dependency graph is 1, so from Theorem 13 . 71 we would obtain 3C 1 ^ C(/). 

Example 6.3. Consider x A y = 2x + 2y defined over F3. It is easy to show that A is commutative and 
idempotent, but not associative. Consider / : F3 — > F|, defined by /(xi,X2) = (2xi + 2x2, xi). It is not 
difficult to show that / has the limit cycle of length 3, {(0, 1), (2,0), (1,2)} (it also has 3 fixed points and 
another limit cycle of length 3). On the other hand, the loop number of its dependency graph is 1, so from 
Theorem [Owe would obtain 3C 1 ^ C(f). 
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7. Discussion 

In this paper we have identified a broad class of discrete dynamical systems with a finite phase space for 
which one can derive strong results about their long-term dynamics in terms of properties of their dependency 
graphs. We classify completely the limit cycles of semilattice networks with strongly connected dependency 
graph and provide polynomial upper and lower bounds in the general case. It is our hope that the formulas 
in this paper are related to general properties of semilattices, which is a subject for future investigation. 
As mentioned in the Introduction, the motivation for this investigation was the need for theoretical tools 
to analyze discrete models in biology. An example of such an application is given in [5J, where it is shown 
that the results about conjunctive Boolean networks can be applied to determining the limiting behavior of 
certain types of epidemiological models. 

The results in this paper apply to certain types of Boolean networks and cellular automata, which in many 
cases have the property that the update functions are of the same type for all nodes. Another model type 
to which the results in this paper apply in some cases is that of so-called logical models, developed by Rene 
Thomas for the purpose of modeling gene regulatory networks. It is shown in [5] that logical models can 
be translated into the framework of polynomial dynamical systems. If the dynamical system arises from a 
semilattice function, then the results of this paper give information about the steady states and limit cycles 
of the model under synchronous update. 
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